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Abstract In the framework of the planar t-J model for cuprates we analyze the development 
of a pseudo gap in the density of states, which at low doping starts to emerge for 
temperatures r < J and persists up to the optimum doping. The analysis is based 
on numerical results for spectral functions obtained with the finite-temperature 
Lanczos method for finite two-dimensional clusters. Numerical results are addi- 
tionally compared with the self consistent Bom approximation (SCBA) results 
for hole-like (photoemission) and electron-like (inverse photoemission) spectra 
at T = 0. The analysis is suggesting that the origin of the pseudo gap is in 
short-range antiferromagnetic (AFM) spin correlations and strong asymmetry 
between the hole and electron spectra in the underdoped regime. 

We analyze also the electron momentum distribution function (HMD). Our 
analytical results for a single hole in an AFM based on the SCBA indicate an 
anomalous momentum dependence of HMD showing "hole pockets" coexisting 
with a signature of an emerging large Fermi surface (FS). The position of the 
incipient FS and the structure of the EMD is determined by the momentum of 
the ground state. The main observation is the coexistence of two apparently 
contradicting FS scenarios. On the one hand, the (5-function like contributions at 
{it/2, n/2) indicate, that for finite doping a pocket-like small FS evolves from 
these points, provided provided that AFM long range order persists. On the other 
hand, the discontinuity which appears at the same momentum is more consistent 
with with infinitesimally short arc (point) of an emerging large FS. 

1. GENERALIZED t-J MODEL 

Spectral properties of underdoped cuprates are of current interest, in particu- 
lar results for the electron spectral functions as obtained with the angle resolved 
photoemission (ARPES) [1, 2, 3]. A remarkable feature of the ARPES data 
is the appearance of a pseudogap already at temperature T* well above the 
superconducting T^. Other quantities, e.g., the uniform susceptibility, the Hall 
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constant and the specific heat also show a pseudogap consistent with energy 
scale T* [4]. At very low temperature T <^ T* the Fermi surface and the 
corresponding electron spectral functions change dramatically with doping of 
planar cuprate systems with holes, where a transition from "small" to "large" 
FS seems to be consistent with ARPES, but is not adequately understood from 
the theoretical point of view. 

There have been also several theoretical investigations of this problem, using 
the exact diagonalization (ED) of small clusters [5, 6, 7, 8], string calculations 
[9], slave-boson theory [10] and the high temperature expansion [11]. While 
a consensus has been reached about the existence of a large Fermi surface in 
the optimum-doped and overdoped materials, in the interpretation of ARPES 
experiments on underdoped cuprates the issue of the debate is (i) why are 
experiments more consistent with the existence of parts of a large FS - Fermi 
arcs or Fermi patches [3, 12] - rather than with a hole pocket type small FS, 
predicted by several theoretical methods based on the existence of AFM long 
range order in cuprates, (ii) how does a partial FS eventually evolve with doping 
into a large closed one. 

The main emphasis of the present study is on the pseudo gap found in ARPES 
and also in some exact diagonalization studies [13, 14, 15]. We employ the 
standard t-J model to which we add a nearest neighbor repulsion term V, 

H = -t (c^W + 

<U>,o- 

+ ^ [JStS^^ + lj{StS- + S-S^) + {V- j)n,nj]. (1.1) 
<ij> ^ ^ 

Here i,j refer to planar sites on a square lattice and cj^ represent projected 
fermion operators forbidding double occupation of sites. Sf are spin operators. 
For convenience we treat the anisotropy 7 as a free parameter, with 7 = in 
the Ising case, and 7 ^ 1 in the Heisenberg model. 

Numerical results presented here were obtained with Lanczos exact diag- 
onalization (ED) technique on small clusters with A'^ = 16 ~ 32 sites, for 
temperature in the range T < J. The method is simple: we take into account 
only the lowest ^ 100 Lanczos states and evaluate the corresponding thermal 
averages. The method is compared with a more elaborate finite-temperature 
Lanczos method (FTLM) [13] and the agreement in here studied T < J regime 
is excellent. It should be noted that for low temperatures the results of the 
diagonalization of small clusters always have to be examined with caution, 
because low energy scale exhibits very strong finite size effects. 

Our analytical approach is based on a spinless fermion - Schwinger boson 
representation of the t-J Hamiltonian [16] and on the SCBA for calculating 
both the Green's function [16, 17, 18] and the corresponding wave function 
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[19, 20]. The method is known to be successful in determining spectral and 
other properties of the quasi particles (QP). In contrast to other methods the 
SCBA is expected to correctly describe the long-wavelength physics, the latter 
being determined by the linear dispersion of spin waves, whereas the short- 
wavelength properties can be studied with various other methods. Here we 
compare the SCBA results with the corresponding ED, as shown further-on. 

In the SCBA fermion operators are decoupled into hole and pseudo spin 
- local boson operators: = h\, Cj^j^ = h\S^ ~ h\ai and Cj^j^ = h\, Ci^i = 
h\S^ ^ h\ai for i belonging to A- and S-sublattice, respectively. The effective 
Hamiltonian emerges 

H = Ar-V2 ^(Mkq4_q/lk< + H.C.) + J2 ^q<«q, (1-2) 
kq q 

where h]^ is the creation operator for a (spinless) hole in a Bloch state. The 
AFM boson operator aq creates an AFM magnon with the energy a;q, and Mkq 
is the fermion-magnon coupling. 

We calculate the Green's function for ahole 0^(00) within the SCBA [16, 17, 
18]. This approximation amounts to the summation of non-crossing diagrams 
to all orders and the corresponding ground state wave function with momentum 
k and energy eu [19, 20, 21] is represented as 

l*k) = 4/'[4 + - + Ar-"/^E^kq,Gk,(a;i)...M£„_^q„x 

qi.--,qn 

X G^Ju;n)hl^a{^...a{^+ ... ]|0). (1.3) 

Here = k— qi— q^, Um = Ck— ^^qi— (^q^ and Zk is the QP spectral 
weight. 

2. PSEUDO GAP IN THE DENSITY OF STATES 

We study here the planar density of states (DOS), defined as M{io) = 
2 / N J2k. — fi), where Ak(i^) is the electron spectral function [14, 13], 
and fj. denotes the chemical potential. First we calculate the DOS with the finite- 
temperature Lanczos method for clusters of iV = 18, 20 sites doped with one 
hole, Nfi = 1. Here we denote with A/^~(a;) the density of states corresponding 
to adding a hole into the system and thus to the photoemission experiments, 
while M^{u!) represents the inverse photoemission (IPES) spectra. 

In Fig. 1.1(a) we present M{uj) for two different J/t = 0.3 and J/t = 0.6 
on a = 18 sites cluster for V = Q. We observe that the pseudo gap scales 
approximately as 2 J. The analysis at elevated temperatures shows that the gap 
slowly fills up and disappears at T ~ J. In Fig. 1.1(b) results for different 
values V are presented. The gap remains robust also in the presence of the 
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Figure 1.1 M{uj) for one hole on N = Id, sites, presented for different J/t, V/t and T/t. (a) 
1/ = 0. (b) J = 0.3t. Broadening of peaks is taken 6/t = 0.04. 

V term, which enhances (V < 0) or suppresses (V > 0) the binding of hole 
pairs. It is thus evident that the effect of V is only of qualitative character and 
not relevant for the existence of the pseudo gap. We therefore believe that this 
analysis suggests that the origin of the pseudo gap is in short-range AFM spin 
correlations rather than in the binding tendency of doped holes. 

AFM spin correlations are correctly taken into account in the SCBA. There- 
fore for the limiting case of low doping, Ch ^ 0, and T ^ we approximate 
A/"" (a;) with 

Ar-(a;) cx^ImGk(-w), (1.4) 

k 

where Gk(w) is the SCBA Green's function for adding one hole (ARPES) to 
an AFM reference system (instead of adding one hole to the state with one 
hole). The corresponding DOS for removing a hole (IPES) from the state with 
one hole, J\f~^{u:) = ^ Xik ^k l"^)' calculated accurately in the SCBA 

as follows. First the spin averaged hole-like spectral function, 

^kM = -^Im^(M/ko|ck,a^^4,al*ko), (1.5) 
■^'1 ^ UJ — ri 

is expressed in terms of holon and magnon operators and the result for M'^ (a;) 
emerges, 

■^^(^) = ^ E(*ko|/il[<^('^ -H)+ ai5{u; - H)a'^^hi\^^,,). (1.6) 

i 

Here |^ko) represents a weakly doped AFM, i.e., it is the ground state (GS) 
wave function of a planar AFM with one hole and the GS wave vector ko. 
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Figure 1.2 M{lo) for iVh = 1 on = 20 sites, with J/t = 0.3, V = Q, T/t = 0.05 
(full line). Dashed heavy line represents the SCBA result on large lattice obtained as a sum 
of M^{io) and M^(uj). The SCBA result is obtained on a AT = 64 x 64 cluster and for 
undoped reference system. Note the "string states" resonances, absent in the finite doping 
Green's function. Af^ (uj) corresponds to IPES. Reference hole concentration is Ch = \/N. 
The SCBA result is normalized to Ch = 1/20. Broadening of peaks is taken 6/t = 0.01. 



The normalization (sum rule) of A^{uj) and M'^i^uj) is discussed in detail in 
Ref. [13]. It should be noted that the normalization of A^{lo) is not trivial and 
is related to the EMD presented in the next section. 

In Fig. 1.2 are shown spectra J\f{oj) obtained with the ED on a = 20 sites 
cluster. We compare these spectra with the DOS within the SCBA. The peaks 
in J\f'^{iL') can well be explained with magnon structure of single hole ground 
state, while peaks in Af~ (uj) wee string states known in the single hole case 
and are for the present study of the pseudo gap and FS not relevant. As seen 
in Fig. 1.2 the total DOS obtained with the SCBA and in particular M'^{u) 
remarkably accurately resemble the ED DOS. 

3. LARGE OR SMALL FERMI SURFACE? 

The electron momentum distribution function nk = ('I' ko I Z^cr Ck cr'^k,crl^ko) 
is the key quantity for resolving the problem of the Fermi surface. Numeri- 
cally rik can be determined by exact diagonalization of the model in small 
clusters, where only a restricted number of momenta k is allowed. The GS 
wave vector due to finite size effects varies with N. Therefore we present here 
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(-71,-71) (0,0) (7t,7t) (7t,-7l) 

Figure 1.3 N{nu, — 1) obtained from ED for various systems A'^ = 16, 18, 20, 32 and J/t = 
0.3. For the SCBA A?^ = 64 x 64, 7 = 0.999 and note that delta-function contributions at 
k = ±ko are not shown. In plotting the curve for iV = 32 data from Ref . [8] were used. 



results obtained with the method of twisted boundary conditions [22], where 
tjji tjji expidjji. Since nt = nk(ko,^) depends both on kg and 9 it 
follows from Peierls construction that nk(ko, 0) = nk+ko(0) ko) for ^ = ko. 
This allows us to study nk for arbitrary k and kg. Furthermore, the finite size 
effects of the results are suppressed if we fix ko for all clusters here studied to 
the symmetry point ko = (f , f ). 

In Fig. 1.3 we present for J/t = 0.3 ED results for clusters with different 
N and 7 = 1. The EMD obeys the sum rule ^Zk n]f_ = N — 1 and, for the 
allowed momenta, the constraint A^(nk — 1) < 1. We show here the quantity 
N{n\^ — 1), which for different N scales towards the same curve. Results are 
presented for some selected directions in the Brillouin zone (BZ) and should 
be averaged over all four possible ground state momenta when discussed, e.g., 
in connection with ARPES data. 

Analytically we study the EMD again in the spinless fermion - Schwinger 
boson representation. The wave function Eq. (1.3) corresponds to the projected 
space of the model Eq. (1.1) and nk = (^ko l^kl^ko) with the projected 
electron number operator nk = Z^a o-'^k,(7- Consistent with the SCBA 
approach, we decouple nk into hole and magnon operators. 
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Figure 1.4 (a) SCBA result for TV = 64 x 64 and J/t = 0.3, 7 = 0.99. (b) Exact diagonal- 
ization result for AT = 32 sites as in Fig. 1.3. 

where ?7j^ = e-*'^ (^»"^:')(l±e"*Q (^«"^j))/2 withQ = {tt.tt). Local a| are 

further expressed with proper magnon operators Oq. In general the expectation 
value nk for a single hole has the following structure [23] 

= 1 - ^^ko((^kko + ^kko+q) + ;^^"k- (1-8) 

Here the second term proportional to (5-functions corresponds to hole pockets. 
Note that 5nk, for the case of a single hole fulfills the sum rule X^k '^'^k = 
Zko — 1 and 5n\f_ < 1. The introduction of 5n\f_ is convenient as it allows 
the comparison of results obtained with different methods and on clusters of 
different size N. 

In Fig. 1.3 we also present the SCBA result. We have also checked the 
convergence of 5n\^ with the number of magnon Unes, n. For J/t > 0.3 we 
find for all k that the contribution of terms n > 3 amounts to less than few 
percent. This is in agreement with the convergence of the norm of the wave 
function, which is even faster [20]. 

In Fig. 1 .4(a) we present this 5nk for the whole BZ. The important ingredient 
of the SCBA is the gapless magnons with linear dispersion and a more complex 
ground state of the planar AFM. Gk(i^) and ek are strongly k-dependent. As 
a consequence nk in general depends both on k and ko. The ground state is 
for the t-J model fourfold degenerate and we again choose ko = (f , f ). In 
Fig. 1 .4(b) is presented also the result of exact diagonalization of a A?" = 32 
sites cluster [8], but generalized to the whole BZ. 

The main conclusion regarding the EMD is the coexistence of two apparently 
contradicting Fermi-surface scenarios in EMD of a single hole in an AFM. (i) 
On one hand, the (5-f unction contributions in Eq. (1.8) seem to indicate that at 
finite doping a delta-function might develop into small Fermi surface, i.e., a 
hole pocket, provided that AFM long range order persists, (ii) A novel feature 



8 



is that also Suk is singular in a particular way, i.e., it shows a discontinuity 
at k = ko with a strong asymmetry with respect to ko. It is therefore more 
consistent with infinitesimally short arc (point) of an emerging large FS. For 
finite doping the discontinuity could possibly extend into such a finite arc (not 
closed) FS. Note that as long-range AFM order is destroyed by doping, hole 
pocket contributions should disappear while the singularity in 5nk could persist. 
The results of the two methods, the SCBA and the ED agree quantitatively at 
all points in the BZ. However, the SCBA result is symmetric around T point in 
the direction k || ko, while small system results show a weak asymmetry for 
k = ±ko, respectively. From our analysis of the SCBA results for — > oo 
and long range AFM spin background it follows that in the thermodynamic 
limit Ch — > nk is symmetric. The asymmetry is in Ref. [8] attributed to 
the opening of the gap in the magnon spectrum at q ~ Q in finite systems. 
Within the SCBA the asymmetry also appears if the HMD is evaluated with ko 
displaced from {^,^)by a. small amount (5ko (not shown here). 

4. SUMMARY AND ANALYTICAL RESULTS 

Full numerical results are captured with a simple analytical expansion which 
gives more insight into the structure of ^4^ (w) and Srik. We simplify the wave 
function, Eq. (1.3), by keeping only the one-magnon contributions and take 
J/t^l and the leading order contributions are then 

Sn^^ = --^koAfkoqGko(eko-'^q)[2^iq + MkoqGko(eko-Wq)] 

~ -8ZlJ^{l + Z^^), q^O, (1.9) 

'^q "^q 

with q = k — ko (or k — ko — Q) and v = t(sin fcoa;) sin koy). 

A surprising observation is that the EMD exhibits for momenta k ~ ko , ko + 
Q a discontinuity ~ ZkoN^^'^ and Sn^^^ oc -(1 + sign qx)/qx- These discon- 
tinuities are consistent with ED results. One can interpret this result as an 
indication of an emerging large Fermi surface at k ~ ikg. The discontinuity 
appears only as points =bko, not lines in the BZ. Note, however, that this re- 
sult is obtained in the extreme low doping limit, i.e., Ch = 1/iV and it is not 
straightforward to generalize it to the finite doping regime. 

A direct reflection of the anomaly in (5nk is also the structure of Ak('^) 
as presented in Fig. 1.5(a), where ^k(^) exhibits a typical asymmetry in 
momentum dependence. This shows the instability towards the large FS. The 
electron part is in this figure approximated with A^{ll!) ~ Z]^^6{e]i — ek,, — to). 

In Fig. 1.5(b) Sn^^^ is presented and the anomaly discussed above is clearly 
seen. 
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Figure 1.5 The tendency of fomation of the large FS as seen in (a) Ak{uj) and (b) 5n]^ ' for 
Zkgt/ J ~ 1. The figures are in arbitrary units and chpped for convenience. 



We conclude by stressing that the origin of the pseudo gap found in cuprates 
seems to be in the short range spin correlations of the reference AFM system, as 
well as in the strong asymmetry between the hole-like and electron-like spectra 
in underdoped systems. From the present SCBA analysis it is clear that the gap 
size is a natural consequence of magnons with dispersion ^ 2 J, and the first 
peak in J\f^{uj) thus corresponds to the "van Hove" high density of magnon 
states. In addition in making contact with ARPES experiments we should note 
that ARPES measures the imaginary part of the electron Green's function. We 
must note that using these experiments in underdoped cuprates rik can be only 
qualitatively discussed since the latter is extracted only from rather restricted 
frequency window below the chemical potential. Nevertheless our results are 
not consistent with a small hole pocket FS (at least only a part of presumable 
closed FS is visible), but rather with partially developed arcs resulting in FS 
which is just a set of disconnected segments at low temperature collapsing to 
the point [1, 3]. The SCBA results for singular Jrik seem to allow for such a 
scenario. It should also be stressed that the SCBA approach is based on the 
AFM long-range order, still we do not expect that finite but longer-range AFM 
correlations would entirely change our conclusions. 
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